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ABSTRACT 

We present a method for correcting near-infrared medium-resolution spectra 
for telluric absorption. The method makes use of a spectrum of an AOV star, ob¬ 
served near in time and close in airmass to the target object, and a high-resolution 
model of Vega, to construct a telluric correction spectrum that is free of stellar 
absorption features. The technique was designed specihcally to perform telluric 
corrections on spectra obtained with SpeX, a 0.8 — 5.5/im medium-resolution 
cross-dispersed spectrograph at the NASA Infrared Telescope Facility, and uses 
the fact that for medium resolutions there exist spectral regions uncontaminated 
by atmospheric absorption lines. However, it is also applicable (in a somewhat 
modihed form) to spectra obtained with other near-infrared spectrographs. An 
IDL-based code that carries out the procedures is available for downloading via 
the World Wide Web. 

Subject headings: techniques: spectroscopic — methods: data analysis — atmo¬ 
spheric effects 


1. Introduction 

Ground-based near-infrared spectroscopy (~ 1—5 pm) has always been hampered by the 
strong and variable absorption features arising from the Earth’s atmosphere. Even within 
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the well-established photometric bands such as J, H, and K, so-called “telluric absorption 
features” are present. The common method used to correct for such absorption features, and 
to retrieve the signal of the object under observation, is to observe a “telluric standard star”, 
near in both time and sky position (airmass) to the object, and subsequently to divide the 
object’s spectrum by the standard star’s spectrum during the hnal phases of the reduction 
process. 

A-type stars are frequently used as such telluric standards, as their spectra contain rela¬ 
tively few metal lines, with strengths that are generally only a few percent of the continuum. 
Furthermore, for A stars with little reddening, the magnitudes at all wavelengths can be de¬ 
duced from a single (usually h^-band) measurement. Finally, the near-infrared continua of A 
stars can be reasonably well approximated by a blackbody with a temperature of ~ 10000 K. 
These latter two properties make A-type stars very useful for determining the absolute flux 
calibration of the observations. However, the strong intrinsic photospheric hydrogen (H) 
absorption features in the spectra of A-type stars present a substantial problem for their 
use as telluric standards, particularly if one is working in the J, ; 2 , Y, or H bands (where 
the numerous higher level lines of the Paschen and Brackett series are located), or if the 
observations are designed to study the H lines in the spectrum of the object of interest. One 
technique commonly adopted in order to preserve the other desirable qualities of A stars 
as telluric and flux standards is to simply interpolate over the H absorption lines in their 
observed spectra. However, since atmospheric absorption features are certainly present at 
the locations of most of these lines, this method can lead to uncertainties in the reality of any 
features subsequently detected in the corrected object’s spectrum at or near the wavelengths 
of the stellar H lines. Furthermore, A-type stars can have rather large rotational velocities 
(upwards of 300 km s“^; e.g., Abt & Morrelli 1995; Royer et ah 2002), which broaden the 
absorption lines. The large number of such broad lines located near the various H series 
limits renders the interpolation method infeasible. 

Maiolino, Rieke, & Rieke (1996) developed a correction method using G2V stars, rather 
than A-type stars, as telluric standards. Their method makes use of the high signal-to-noise 
(S/N), high-resolution spectrum of the Sun, obtained using a differential airmass technique at 
the McMath Telescope at Kitt Peak National Observatory by Livingston & Wallace (1991; 
see also Wallace, Hinkle, & Livingston 1993), to generate a telluric correction spectrum 
which is then applied to the spectrum of the object of interest. The correction spectrum is 
derived from the ratio of the solar spectrum (shifted to the appropriate radial velocity and 
convolved to the resolution of the observations) to the observed G2V star spectrum. This 
ratio should, in principle, remove any intrinsic stellar features from the telluric correction 
spectrum. Although Maiolino et al. (1996) demonstrated that this method can be reasonably 
successful in performing telluric corrections on moderate S/N and resolution data, there are 
several practical problems with it, especially when applying it to higher quality and resolution 
spectra. Firstly, the number of reasonably bright G2V stars available for use as telluric 
standards near any given object is less than the number of A-type stars, simply because 
G stars are inherently less luminous than A stars and catalogues of stars with well-dehned 
spectral classifications are based on magnitude-limited samples. This is demonstrated in 
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Table 1, in which we have compiled the numbers of stars listed with AOV and G2V spectral 
types in the SIMBAD database as a function of V magnitude. The intrinsic (V — K) color 
of a G2V star is about 1.5 mag; therefore, even in the K band, bright AOV stars outnumber 
G2V stars in the database. Secondly, G-type stars have a large number of (relatively weak) 
metal lines throughout their spectra. Hence the accuracy of the resulting corrections depends 
on how well the observed spectrum of the G star matches that of the Sun; the higher the 
S/N and resolution of the data, the closer the spectral match must be for a desired level of 
accuracy. As the strengths of the metal lines vary considerably with spectral type and stellar 
parameters, any mismatch, such as might be caused by spectral classification inaccuracies, 
metal abundance differences, anomalous line strengths, etc., will introduce spurious features 
into the telluric correction spectrum. Thirdly, the reference solar spectrum was obtained 
at a relatively low site and has large wavelength gaps where the intrinsic spectrum could 
not be recovered using the differential airmass method. Because the depths of atmospheric 
absorption features vary with elevation, this solar spectrum would not provide a useful 
correction template for spectra, such as those obtained at a high site such as Manna Kea, 
which have recoverable flux in the regions of strong absorption between the photometric 
passbands. 

An alternative method, used for example by Hanson, Gonti, & Rieke (1996), combines 
the A-star technique with that of Maiolino et ah (1996) and involves the following steps: 
observe both A-type stars and G2V stars, normalize the spectra of both stars, “correct” 
the normalized G2V star spectrum by dividing it by the solar spectrum, replace the regions 
centered on the H lines in the normalized A star spectrum by the corrected G star spectrum, 
and finally divide the result into the object’s spectrum. This is equivalent to dividing the ob¬ 
ject spectrum by that of the A star, and multiplying the result by a line correction spectrum 
generated from the G2V star and the solar spectrum. Since the solar spectrum is used only 
in the regions of the H lines, the resulting correction template should be more accurate than 
one based on the G2V star/solar spectrum ratio alone. Because the normalization process 
removes the signature of the instrumental throughput from the observed stellar spectra, the 
object must be flux-calibrated separately, after the correction process, usually by multiply¬ 
ing the normalized, corrected spectrum by a blackbody with a temperature appropriate for 
an A star. While this method might seem attractive in principle, in practice it requires 
observations of two sets of standards (a suitable one of which may be difficult to hud, as 
explained above) for each object of interest (ideally) and involves a somewhat tedious and 
complicated reduction procedure, particularly if the wavelength range of interest includes 
several absorption lines of the Pa or Br series. 

Glearly, the ideal telluric standard would be one whose spectrum is completely fea¬ 
tureless and precisely known. In addition, such objects should be reasonably bright, fairly 
common, and widely distributed on the sky, so that one can be found close to any particular 
object of interest. While few objects satisfying the first criterion are known in nature, there 
are a limited number of objects for which precise and accurate theoretical models are avail¬ 
able. Vega is one such object. In this paper we describe a new method of generating a telluric 
correction spectrum using observations of AOV stars and a high-resolution model of Vega’s 
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intrinsic spectrum. Through a process somewhat similar to that described by Maiolino et 
ah (1996), this technique preserves the desirable qualities of A-type stars, while avoiding 
the problems with the intrinsic stellar absorption lines, and does not require observations of 
G2V stars. 

The method outlined below was developed for the reduction of spectra obtained with 
SpeX at the NASA Infrared Telescope Facility (IRTF) on Manna Kea. SpeX is a medium 
resolution (A/AA ~ 200 — 2500) cross-dispersed spectrograph equipped with a 1024 x 1024 
InSb array which provides simultaneous wavelength coverage over 0.8 —2.4 /rm in one grating 
setting and 2.4 — 5.5 /rm in another grating setting (see Rayner et ah 1998, 2002 for a full 
description of SpeX and its various scientihc observing modes). Because the entire Paschen 
series and most of the Brackett series of H absorption lines are simultaneously detected in 
a spectrum of any AOV star obtained with SpeX, it was important to devise a technique 
that could be applied easily and rapidly to accurately remove all of these features from 
a telluric correction spectrum generated from the AOV star spectrum. With high quality 
spectra of an AOV star, a telluric correction spectrum accurate to better than about 2 
percent across the H lines can be quickly constructed with this technique. Although designed 
specihcally for SpeX data, the method can also be used to correct spectra obtained with 
most medium resolution (A/AA < 50000) near-infrared spectrographs. A set of IDL-based 
routines (SpextooLextension, v. 2.0), specihcally designed for application to SpeX data, 
implements the telluric correction method described below and comprises an extension to 
Spextool (v. 2.3), a package of routines written for the basic reduction of spectra obtained 
with SpeX (Cushing, Vacca, & Rayner 2002). Both sets of programs are available from the 
SpeX webpage at the IRTF website (http://irtfweb.ifa.hawaii.edu/Facility/spex). 
Finally, a third set of IDL routines, constituting a general version (i.e., not specihc to SpeX 
data) of the telluric correction code, is also available from this webpage. 


2. Description of Method 

We assume that the observed spectrum 0(A) of any object is given by 

0(A) = |/(A).T(A)]»P(A).Q(A), (1) 

where /(A) is the intrinsic spectrum of the object, T(A) is the atmospheric (telluric) absorp¬ 
tion spectrum, P{X) is the instrumental prohle, Q{X) is the instrumental throughput, and 
* denotes a convolution. Because the convolution represents a smoothing of the intrinsic 
stellar and telluric spectra, the above equation can be re-written as 

0(A) = |/(A) » P(A)|. |T(A) » P(A)|. Q(A), (2) 

or 

0(A)= W(A)-R(A), (3) 


where 


^instr (A)=J(A)*P(A) 


(4) 
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is the spectrum of the object observed above the atmosphere with an instrument with perfect 
throughput, and 

fi(A) = |r(A)»P(A)].Q(A) (5) 

is the telluric absorption spectrum convolved with the instrument prohle and multiplied by 
the instrumental throughput. The latter represents the overall throughput, or response, of 
the atmosphere+instrument system. We assume that Q(A) is a smooth, featureless function 
and therefore any observed high frequency variations in 0(A) must be due to either 1(A) or 
T(A). 

The immediate objective is to determine Unstr of a specihc object from 0(A)/i?(A). 
Conversely, if /(A) for any particular source is known a priori, and P{X) of the instrument 
is known or can be determined, then the system throughput R{X) can be determined from 
Eqs. (3) and (4). This factor can then be used to determine linstr(^) of a target object from 
the observed spectrum, provided that the telluric absorption spectrum T(A) does not vary 
substantially between the two sets of observations. This forms the basis of both our method 
and that of Maiolino et ah (1996): observations of a “telluric standard”, whose intrinsic 
spectrum is known, to derive a system throughput curve which can then be applied to a 
target object spectrum observed nearby in the sky and close in time. Since T(A) varies with 
airmass, and on time scales of the order of several to tens of minutes (depending on the 
atmospheric conditions), it is clearly best to observe a tellnric standard at an airmass as 
close as possible to that of the target object, and within a few minntes of the observations 
of the target. 

Because the metal lines in AOV spectra are relatively infrequent and weak (on the order 
of a few percent of the continnnm level in the near-infrared at the resolntion of SpeX; see, for 
example. Fig. 1), the intrinsic near-infrared spectra of snch objects consist essentially of H 
absorption lines and a smooth, featnreless (blackbody) continnnm. Thus, metal abundance 
variations do not strongly affect the appearance of the observed spectra and, except for 
possible differences in the widths and depths of some of the H lines (dne to different rotational 
velocities or slightly different valnes of the surface gravity), the spectra of all AOV stars can 
be considered nearly identical to one another. It is a fortnnate circumstance that Vega, 
the object that defines the zeropoint of the magnitnde system, is also a bright AOV star 
whose spectrnm has been accnrately modeled. Assnming the model provides an accnrate 
representation of the intrinsic spectrnm of Vega, and that Vega is a good representative (an 
archetypal example) of the class of AOV stars, we can use the theoretical model to determine 
Anstr(A) of any AOV star. Hence, if each observation of a given target object is preceded 
and/or followed by an observation of an AOV star, and if the model spectrnm of Vega can be 
scaled and reddened to match the near-infrared magnitndes of the observed AOV star, and 
modified to acconnt for the differences in line strengths, radial and rotational velocities, and 
spectral resolntion, we will be in the happy sitnation of having at hand both the observed 
spectrum, 0(A), and intrinsic spectra, /(A) and linstr(A), of the AOV star. The system 
response fnnction, R{X), can then be determined immediately and used to correct the target 
object’s spectrum. The difficulty then lies “only” in properly modifying the model spectrum 
of Vega. 
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Modifying the model Vega spectrum involves shifting it to the radial velocity of the 
observed AOV star, scaling and reddening it to reproduce the known magnitudes of the 
observed AOV star, convolving it with a function (kernel) that broadens the lines to the 
observed widths and smooths it to the observed resolution, resampling the result to the 
observed wavelength scale, and altering the depths of various H lines to match those of the 
observed AOV star. The convolution kernel is given by the instrument prohle (IP), P’(A), 
and a (rotational) broadening function 0(A) (see below). There are at least three ways 
of constructing T’(A): (1) assume a theoretical prohle whose properties are dehned by the 
instrument parameters (e.g., slit width); (2) construct an empirical prohle using observations 
of unresolved arc lines (similar to those used to perform the wavelength calibration of the 
spectra); and (3) construct a prohle using the observations of the AOV star itself. We 
discuss the last method in detail here, and describe below a set of routines that implement 
this method for use in the reduction of data obtained with SpeX. (The routines also allow 
users to adopt empirical or theoretical instrument prohles.) 

We begin by selecting a wavelength region centered on an intrinsic absorption line in the 
observed and model AOV spectra. The observed spectrum is normalized in the region of the 
line by htting a low order polynomial or a spline to the continuum; a good representation 
of the observed continuum can be obtained once bad pixels and the regions around any 
strong absorption features are excluded from the htting process. The normalization removes 
the signature of the instrumental throughput as well as the shape of the stellar continuum. 
Once the observed stellar spectrum has been normalized, the wavelength region around the 
absorption feature can be cross-correlated with that from the normalized model of Vega, 
using a technique similar to that described by Tonry & Davis (1979), to determine the 
observed radial velocity of the AOV star. The model is then shifted by the estimated radial 
velocity, and scaled and reddened to match the observed magnitudes of the AOV star. 

We now wish to construct a convolution kernel which, when applied to the model spec¬ 
trum of Vega, will produce a close approximation of the intrinsic spectrum of the observed 
AOV star, as measured by an instrument with perfect throughput. That is, we assume that 
the intrinsic spectrum of the AOV star can be given by 

IaovW = Ivegai^) * ©(A) (6) 

where IvegaW is the scaled and shifted model spectrum of Vega and 0(A) is a function 
that accounts for the differences in rotational velocity and any other line broadening factors 
between the Vega model and the observed AOV star. The observed AOV spectrum Oaov(A) 
is then given by 


Oaov(A) = {[IvegaW * ©(A)] ■ T(A)} * P(A) ■ g(A), (7) 

or, 

OaOvW = [IvegaW * K{\)] ■ P(A), (8) 


where 


K{\) = 0(A) *P(A), 


(9) 
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and 

Iinstr,AOV Ivega{X)*K{X). (10) 

For the moment, let us simply assume we can find an intrinsic stellar absorption line in the 
observed AOV spectrum across which both the atmospheric absorption and the instrumental 
throughput are either constant or only smoothly varying. Using the normalized observed 
and model spectra in the vicinity of the line, we then have, 

OW(A)=/ve,a(A)*i^(A), (11) 

where primes denote the normalized spectra (divided by the respective continua). The 
convolution kernel K{X) can then be determined by selecting a small region around the 
observed absorption line and performing a deconvolution. 

Once K{X) has been determined, it is convolved with the shifted, scaled, and reddened 
Vega model spectrum and R{X) is determined from equation (8). Under the assumptions 
that all AOV stars are identical to Vega and that the model of Vega provides a highly ac¬ 
curate representation of Vega’s intrinsic spectrum, the convolution and subsequent division, 
OAovW/[Ivega{X) * A'(A)] should yield a response spectrum in which the intrinsic stellar 
absorption lines are absent. In practice, we find that some small residual features at the 
locations of the H lines are discernible in this response spectrum due to differences in the 
strengths of these lines between the model and the observed AOV spectra. Some of these 
differences may in fact be real, due to different values of logg between the standard star 
and the Vega or misclassification of the observed standard star and hence a slight spectral 
mismatch between the standard and Vega. However, the model H line strengths (equivalent 
widths) can be easily adjusted to minimize, or eliminate, these residuals. Furthermore, since 
any metal lines are fairly weak and do not vary strongly among late B-type and early A- 
type stars, the effects of spectral mismatches and misclassifications on the resulting telluric 
correction spectrum can be minimized to a large extent by simply modifying the H line 
strengths in the model. For the resolving power and typical S/N levels achieved with SpeX, 
minor misclassifications should have a negligible effect on the final corrected spectra. 

The intrinsic spectrum of the target object linstrW is then determined from equation 
(3), re-written as 

RnstriX) = 0(A) ■ [IvegaW * K{X)/OaovW]. (12) 

We shall refer to the quantity [Ivega{X) *K{X)]/Oaw{.^) = 1/-R(A) as the “telluric correction 
spectrum”. One advantage of this method is that the resulting spectrum Iinstr{X) is auto¬ 
matically placed on the proper flux scale, if the slit losses during the two sets of observations 
(standard star and target object) were the same. 

In summary, the method of constructing a telluric correction spectrum consists of an 
observation of an AOV star and several steps designed to modify the model spectrum of Vega 
appropriately: 

1. Normalization of the observed AOV star spectrum in the vicinity of a suitable absorp¬ 
tion feature (as defined below); 



2 . Determination of the radial velocity shift of the AOV star; 

3. Shifting the Vega model spectrum to the radial velocity of the AOV star; 

4. Scaling and reddening the Vega model spectrum to match the observed magnitudes of 
the AOV star; 

5. Construction of a convolution kernel from a small region around an absorption feature 
in the normalized observed AOV and model Vega spectra; 

6 . Convolution of the kernel with the shifted, scaled, and reddened model of Vega; 

7. Scaling the equivalent widths of the various H lines to match those of the observed 
AOV star. 

Finally, the convolved model is divided by the observed AOV spectrum and the resulting 
telluric correction spectrum is multiplied by the observed target spectrum. 


3. Practical Considerations 

For the method described above to work properly, the convolution kernel K{X) must be 
determined accurately. If P{\) and 0(A) are known in advance, K{\) can be constructed 
using equation 9. If K{X) is to be generated from the observations themselves, it is necessary 
to hnd a region of the observed AOV spectrum which contains a strong intrinsic stellar 
absorption line (so that the radial and rotational velocities and IP can be determined), 
but is free from atmospheric absorption or instrumental features. The latter requirement 
arises from our assumption that T(A) and Q{X) are either constant or vary only slowly 
and smoothly across the width of the line. Fortunately, for observations with SpeX, the 
atmosphere has cooperated in providing such a wavelength region and stellar atmospheres 
have cooperated in providing a suitable absorption line. Fig. 1 presents a model of the 
atmospheric transmission between 0.9 and 1.2 /rm, computed with the ATRAN program 
(Lord 1992) for an airmass of 1.2 and standard conditions on Mauna Kea and smoothed to a 
resolving power of 2000. The normalized model spectrum of Vega obtained from R. Kurucz 
(http://kurucz.harvard.edu/stars.html), also smoothed to resolving power of 2000 , is 
overplotted. It can be seen that the Pa S (1.00494 fim) absorption line in the Vega spectrum 
is located in a region where, at this resolution, atmospheric absorption is almost completely 
negligible. Therefore, this line can be used to determine the convolution kernel K{X) in the 
procedure described above. 

In addition to the lack of telluric contamination, the Pa S line also lies in a wavelength 
region where the throughput Q{X) of SpeX is relatively high and does not change dramatically 
across the width of the line. The throughput around the Pa 7 (1.0938 fim) line, for example, 
drops quickly across the line (due to the grating blaze function) and atmospheric absorption 



affects the shape of the long wavelength wing, making the determination of the continuum 
near the wings difficult. Nevertheless, although the Pa S line is ideal for our purposes, 
we have successfully generated telluric correction spectra using the method described in 
Section 2 with kernels built from other H lines in the spectra of AOV stars (e.g., Br 7 ). For 
spectrographs that do not contain cross-dispersers and/or yield only a limited wavelength 
range, the observed AOV spectra may not contain a line that is suitable for constructing K{X). 
In these cases, the convolution kernel can be generated from the IP and some estimate of the 
rotational velocity of the AOV star (Eq. 9). At resolving powers similar to those of SpeX, 
the latter generally has a negligible effect on the resulting telluric correction spectrum. 

The model of Vega we use for our telluric correction procedure has a resolving power 
(A/A A) of 500000, binned by 5 (giving an effective resolving power of 100000), and has been 
scaled to match the observed Vega flux at 5556 A given by Megessier (1995). (A higher 
resolution model is also available from the Kurucz web page; however the array lengths 
necessary to store these data and the subsequent intermediate products become inordinately 
large and unwieldy.) The properties of the model are as follows: Tes = 9550 K, loggf = 3.950, 
Urot = 25 km/s, Uturb = 2 km/s, and a metal abundance somewhat below solar. The model 
also has an associated continuum spectrum; using the full spectrum and the continuum, we 
generated a normalized line spectrum. The continuum spectrum contains sharp jumps due 
to the abrupt increases/decreases in the opacity at the limits of the H (Pa and Br) series. 
The strengths of these jumps must be modihed along with the line strengths. To account 
for this, we have normalized the continuum spectrum by dividing it by a spline htted to 
the spectrum far from the jumps. The £t passes approximately through the flux midpoints 
between the continuum levels on either side of the jumps. Scaling this normalized continuum 
by an appropriate value then has the effect of scaling the heights of the jumps. 

Note that because the Vega model spectrum has such a large resolving power, the 
method we have outlined should work successfully for spectral data with considerably larger 
resolving powers than that of SpeX (A/AA = 200 — 2500). In principle, use of the routines 
described below (Sec. 4) is limited to data with A/AA < 50000, although data with higher 
resolving powers could be corrected if a model with higher resolving power were incorporated 
into the code. We have used a more general version of the code described below to successfully 
perform telluric corrections on data with resolving power of 5000 and we routinely use 
the code, with an empirical instrument profile T’(A), to successfully correct low resolution 
(A/AA ~ 200) prism data obtained with SpeX. At low resolving powers, the possibility 
of finding isolated H lines that are uncontaminated by atmospheric absorption is greatly 
diminished and precludes the construction of the convolution kernel K{X) from the observed 
data. At A/AA ~ 200, for example, even the Pa 6 line is affected by telluric absorption in 
the wings. 

The calculations described above (Sec. 2 ) are easiest to carry out in velocity coordinates, 
which is appropriate for a spectrum with constant spectral resolving power. In addition, the 
Vega model has a constant resolving power. However, in most spectrograph designs, the 
resolving power varies across the spectral range. SpeX is no exception, and the resolving 



power varies linearly by about 20% across each order (Rayner et al. 2002). Fortunately, the 
spectral resolution AA of SpeX is constant to a high degree across each order, although the 
specihc value varies from order to order. However, in pixel coordinates it is the same value 
for all orders, given approximately by the projected slit width. This implies that, in pixel 
coordinates, a kernel K{X) built from a line in one order can be applied to all other orders. 
Therefore, the code described below carries out the deconvolution in wavelength coordinates 
in order to construct the kernel, converts the resulting K{X) to pixel coordinates, and then 
applies this kernel to each order. This saves the user from the tedious task of building kernels 
specihc to each order. 

Deconvolution procedures are notoriously sensitive to noise. To minimize the effects 
of noise on the deconvolution needed to construct K{X), we multiply the ratio of the 
Fourier transforms of 0'j^Qy{X) and I'y^gA^) ^y a tapered window function of the form 
1/(1 + {v/VmaxY^)) where Vmax D a reference frequency given by a multiple of the width 
of a Gaussian £t to the Fourier transforms. This has the effect of suppressing high frequency 
noise. 

Clearly, as the response spectrum R{X) will be divided into the observed target spectrum, 
it is also important to acquire high S/N values in the observed spectrum of the AOV star in 
order to achieve accurate and reliable telluric corrections. In general, we have found that a 
S/N of > 100 is needed in order to obtain telluric-corrected spectra that are limited by the 
hat-helding errors and the S/N of the observed target spectrum. Of course, the accuracy 
of the dual corrected spectra also depends on the stability of the atmospheric conditions, 
the proximity (in both time and airmass) of the AOV observations to those of the target 
object, and the degree to which the AOV spectrum matches the model spectrum of Vega. 
However, as we have argued above, small spectral mismatches should have a minor effect 
on the generated telluric correction spectrum. Fortunately, S/N values of 100 are rather 
easy to achieve at resolving powers of < 2500 on 3-m class telescopes. Using SIMBAD we 
have compiled a list of AOV stars that are suitable for carrying out the method described 
above. The SpeX Webpage includes a web form that will generate a listing of the known 
AOV stars within a specihed angular distance and/or airmass difference from any target 
object. Our experience using SpeX on Manna Kea indicates that on particularly dry and 
stable nights, accurate corrections can be achieved when the airmass differences and time 
intervals between the target and standard observations are as large as 0.1 and 30 to 60 min, 
respectively. However, for typical conditions, good telluric corrections require closer matches 
in airmass and shorter time separations. 

Accurate hux calibrations require knowledge of the near-infrared huxes or magnitudes 
of a set of AOV stars. Unfortunately an accurate and comprehensive set of huxes and/or 
magnitude measurements for AOV stars at these wavelengths is currently not available. 
Therefore, we have chosen to adopt the optical {B and V) magnitudes given by SIMBAD 
and the reddening curve given by Rieke & Lebofsky (1985; see also Rieke, Rieke, & Paul 
1989) to estimate the hux as a function of wavelength for a given AOV star. In order to 
achieve a S/N of about 100 without inordinately long exposure times, the telluric AOV star 
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must be reasonably bright; for SpeX observations, AOV stars with V magnitudes of ~ 5 — 8 
have been found to work best. SIMBAD provides highly accurate magnitudes for a large 
number of AOV stars in this brightness range. In addition, the reddening of such objects is 
necessarily small, even in the visible. Hence the precise reddening curve adopted has a minor 
effect on the results. However, we point out that accurate spectroscopic flux calibrations at 
near-infrared wavelengths would clearly beneht from a dedicated program of careful J, II, 
and K observations of a set of AOV stars. Such a program could also serve as the basis for 
establishing a set of near-infrared “spectroscopic flux standards” similar to those available 
in the optical regime. 

Proper flux calibration also requires that the seeing and other slit losses do not change 
substantially between the observations of the target and the AOV star. If the apertures 
used for the spectral extraction are wide, moderate changes in seeing should not affect the 
resulting flux calibrations greatly. In practice we hnd that for typical seeing conditions on 
Manna Kea at the IRTF and extraction apertures of ~ 2" width, the derived fluxes agree 
with those estimated from broad-band magnitudes to within a few percent. 


4. Implementation: xtellcor 

The method outlined above (Sec. 2) has been implemented in a graphical IDL-based code 
called xtellcor. Here we describe the version of xtellcor specihcally designed to operate on 
SpeX data. A more general version of xtellcor (for non-SpeX data) is similar, but does not 
construct a convolution kernel K{\) from the AOV star; rather the user must input values 
of the parameters describing the IP, H(A). Both versions of the code are available from the 
SpeX Website. Help hies are included in the distribution packages to guide the user through 
the various steps. 

Starting xtellcor brings up an IDL widget containing helds in which the user can enter 
the information necessary for the code to run (data directory and hie names for the observed 
AOV and target object spectra; see Fig. 2). Since near-infrared magnitudes of a large set of 
AOV stars are not available, the user must give the B and V mags of the observed AOV star. 
These magnitudes are used to redden and scale the model spectrum of Vega in order to set 
the hux calibration properly. 

At this stage the user can choose to build a kernel from the observed AOV star, or accept 
an empirical IP constructed from the prohles of arc lines in SpeX wavelength calibration 
images. The latter is particularly useful in cases where the IP cannot be constructed from 
the observed AOV star (e.g., in the long-wavelength cross-dispersed modes of SpeX which 
cover the 1.9-5.5 jum range). The arc line prohles have been htted with an empirical function 
of the form 

P(x) = C ■ [erf((a; -f a)/h) — erf((a; — a)/h)], (13) 

where C* is a normalization constant and a and b determine the width of the prohle in pixel 
{x) coordinates. This functional form provides an excellent ht to the observed prohles of 
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the arc lines over all the orders of SpeX and for slit widths of 0.3", 0.5", 0.8", and 1.6". By 
adopting a £t, rather than the observed line profiles themselves, errors in the IP resulting 
from undersampling, low S/N, etc., are minimized. The best fit values of a and b for each 
slit width are stored in a separate data hie. 

If the user chooses to build the kernel from the AOV star, a second window is generated 
containing the observed spectrum in the requested order. The Pa S line is located in order 
6 of the 0.8-2.4 /rm cross-dispersed mode of SpeX. The user can then interactively select 
continuum regions to be fit with a polynomial of a desired degree. Once an adequate fit 
to the continuum in the vicinity of the spectral line is obtained, the normalized spectrum 
(observed divided by the continuum fit) is displayed (Fig. 3). The user then selects the region 
of the normalized spectrum to use in the deconvolution and the code proceeds to scale both 
the model line spectrum and the normalized model continuum spectrum by the ratio of the 
observed equivalent width to the model equivalent width of the selected line (Fig. 3). The 
relative velocity shift between the observed spectrum and the model is determined and the 
kernel is generated via deconvolution after the model has been shifted. Our experience has 
shown that with careful choice of the normalization and deconvolution regions, a kernel can 
be routinely constructed that yields maximum deviations of the convolved model from the 
observed line profile of less than 2% over the selected wavelength range. We typically achieve 
maximum deviations of less than 1.5% and RMS deviations of less than 0.75% across the Pa 
6 line for data with A/AA = 2000 and S/N > 100. 

The observed spectrum of the AOV star is divided by the line model after the latter has 
been convolved with the kernel (either the semi-empirical IP, or one built from the observed 
AOV star) and the result is displayed for each order (Figs. 4 and 5). This is the first version 
of the response spectrum. At this stage the user can adjust the velocity shift of the model 
and modify the strengths of the various H lines to better match those of the observed AOV 
star. Mismatches in line strengths and velocities are easily seen in the telluric correction 
spectrum as “absorption” or “emission” features at the location of the H lines. The user can 
adjust the individual H line strengths of the model by varying the scale factors applied to the 
model line equivalent widths until residual features are minimized in the response spectrum 
in the regions around the lines (see Figs. 4 and 5). This can be done graphically using the 
mouse until residual features are minimized in the regions around the lines. A representative 
spectrum of the atmospheric transmission is displayed to guide the user in this task. (As 
an additional aid to the user, this atmospheric transmission curve can be multiplied by the 
estimated throughput curve of SpeX.) The code also includes a procedure to determine the 
scale factor for a single isolated line automatically. This procedure requires two points on 
either side of a line, fits a Gaussian plus a continuum to the response spectrum and adjusts 
the strength of the model line until the height of the Gaussian is minimized. Most of the 
H lines can be easily removed from the response/telluric correction spectrum at this stage. 
However, we have found that the Pa 12, Pa 13 and Pa 14 lines (and occasionally the Pa 11 
and Pa 15 lines) almost always leave residual features in the telluric correction spectrum that 
cannot be eliminated, no matter how the line strengths are scaled. This indicates that the 
Vega model we are using does not provide a good match to the prohles of these lines in the 
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spectra of real AOV stars. Once an improved Vega model is available, we will incorporate 
it into the code. Fortunately, even in the vicinity of these discrepant lines, the residual 
featuress are generally less than a few percent and therefore do not result in large errors in 
the hnal spectrum. 

The resulting set of line strength scale factors are joined with a spline £t to generate 
a line scaling spectrum which is multiplied by the line model. The normalized continuum 
(divided by the spline fit) is also multiplied by this line scaling curve. This has the effect 
of modifying the jumps seen in the continuum model. Since the highest order lines rarely 
require modification of their strengths, the jumps are normally scaled by the equivalent width 
ratio determined from the line used to build the kernel. The scaled line model is then used 
to regenerate the response spectrum. When the user is satished with the modihcations to 
the model, the hnal telluric correction spectrum is generated. This can then be shifted in 
pixels in order to align the telluric absorption features with those seen in the target object 
spectrum. This procedure should minimize any wavelength shifts due to hexure between the 
observations of the target object and the standard star. The code automatically computes 
the best pixel shift using a region of the spectra selected by the user (Fig. 6). The shifts 
are usually less than 0.5 pixels for standard stars observed within 10 degrees of the target 
object. The shift value is then applied to all the orders of the telluric correction spectrum. 
Finally, the observed target object spectrum is multiplied by the shifted telluric correction 
spectrum, order by order, and the hnal corrected object spectrum is generated (Fig. 7). The 
user has a choice of units for the hnal output spectrum. The user can also choose to save both 
the telluric correction spectrum itself, as well as the modified Vega spectrum. The former 
is useful if a single AOV star has been observed as a telluric standard for multiple target 
objects; a separate routine allows this previously constructed telluric correction spectrum to 
be read in and used to correct the observed spectra of other objects. 


5. Application and Examples 

The response spectra generated by xtellcor for the various short-wavelength cross- 
dispersed orders of SpeX are shown in Fig. 8(a-f) and Fig. 1 of Cushing et ah (2002a). As 
can be seen from these figures, the features seen in the response spectra match extremely 
well with those seen in the theoretical atmospheric spectrum, and the stellar H lines seen 
in the observed AOV spectra are removed to a very high accuracy. Small residuals can be 
seen in Fig. 8(e) at the locations of the Pa 11 - Pa 15 lines. Aside from these features, 
however, the accuracy of the telluric corrections is limited only by the S/N of the input AOV 
spectrum. The latter is determined by the exposure time, hat-fielding errors, instrument 
throughput, and the atmospheric transmission itself. (Clearly, the derived telluric correction 
spectrum will be highly uncertain in those wavelength regions where the atmosphere is so 
opaque that the observed spectra have S/N close to zero.) Of course, as explained above, 
the accuracy of the hnal corrected target spectra also depends on the diherences in time and 
airmass between the observations of the AOV star and the target object, as well as the S/N 
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of the target spectrum. 

Examples of typical output spectra generated by xtellcor can be seen in Fig. 7 (for HD 
223385, A3Ia), Fig. 9 (for HD 206267, 06.5V((f)), and Fig. 10 (for HD 16139, G8Ha), as well 
as in the papers by Cushing et al. (2002a) and Rayner et al. (2002). The differences in airmass 
and time of observations between the objects shown in these hgures and their telluric AOV 
standards were < 0.05 and < 20 min, respectively, in all cases. However, the atmospheric 
conditions (humidity, seeing, transparency) were variable during the observations. The S/N 
values of the AOV spectra were always > 100. As the hgures demonstrate, the method enables 
us to recover the target spectrum even in regions of relatively low atmospheric transmission. 
We note that the equivalent widths of the intrinsic stellar absorption and emission features 
seen in the spectra in Figs. 7 and 9 agree well with those measured for these spectral classes 
by Hanson et al. (1996) in the K band. 


6. Summary 

We have described a technique for correcting observed near-infrared medinm-resolntion 
spectra for telluric absorption. Designed specihcally for data obtained with SpeX at the 
IRTF, the method makes use of a high S/N (> 100) spectrum of an AOV star observed near 
in time and close in airmass to the object of interest, and uses a theoretical model of Vega 
to remove the intrinsic H absorption featnres seen in the stellar spectrum. The method can 
regnlarly generate tellnric correction spectra accnrate to better than 2% in the vicinity of 
the intrinsic stellar H lines and significantly better than that at other wavelengths. 

We have demonstrated the application of this method on spectra obtained with SpeX. A 
set of graphical, easy-to-use, IDL-based routines (part of a general SpeX spectral redaction 
package called SpextooLextension v. 1.1), implements the procedures described in this paper 
and can be found on the SpeX webpage at the IRTF website 

(http: //irtfweb. if a.hawaii . edu/Facility/spex/). A general version of the code, usefnl 
for performing tellnric corrections on spectra obtained with other near-infrared spectrographs 
with resolving powers of A/A A < 50000 is also available. In addition, a form that generates 
a list of the AOV stars listed in the SIMBAD database near a set of specihed coordinates is 
also posted on the website. 


We thank R. Knrncz for providing the high resolntion model of Vega and J. Tonry for 
suggestions regarding the form of the instrnment prohles. We also thank the SpeX team 
for their expert constrnction of an invalnable instrnment and the staff at the IRTF for their 
assistance at the telescope. This research has made use of the SIMBAD database, operated 
at CDS, Strasbourg, France. 
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Table 1. Cumulative Number of AOV and G2V Stars listed in the SIMBAD Database, as 
of July 2000, as a function of V band magnitude. 


V (mag) 

A(AOV) 

A(G2V) 

< 7 

364 

50 

< 8 

740 

171 

< 9 

1687 

459 

< 10 

3633 

993 

< 11 

4557 

1165 



Atmospheric Transmission 


(Lunnupuoo ‘D50/\)^j/(D60/\)^j 

Fig. 1.— The normalized model spectrum of Vega, smoothed to a resolving power of 2000 
(thick line, left hand axis). The theoretical atmospheric transmission around the Pa 6 line 
(1.00494 fim) at airmass 1.2 for typical conditions on Manna Kea, smoothed to a resolving 

















18 



Fig. 2.— Control panel widget for xtellcor. 
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Fig. 3.— Determination of the convolution kernel from the Pa 6 line in the normalized 
spectrum of an AOV star. The vertical dashed lines denote the region of the spectrum used 
in the deconvolution to generate the kernel. Note that the residuals (seen in the bottom 
panel) have an RMS deviation of much less than 1 % . 
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Fig. 4.— Determination of the H line scale factors for the Vega model in order 4 of SpeX. 
The upper panel shows the scale factors for each of the H absorption lines, while the bottom 
shows the response spectrum. The yellow curve is the theoretical atmospheric transmission 
for an airmass of 1.2 and standard conditions on Mauna Kea. 
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Fig. 5.— Same as Fig. 4 for order 6 of SpeX. Note the small residual features in the response 
spectrum at the locations of Pa 12 and Pa 13. 
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Fig. 6.— Determination of the shift between the telluric correction spectrum and the ob¬ 
served target object spectrum. The vertical dashed lines denote the region used to determine 
the pixel shift. The solid white line is the observed object spectrum, while the green line is 
the system response (the inverse of the telluric correction spectrum). The lower panel shows 
the telluric-corrected spectrum for this order. 
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Fig. 7.— The final telluric-corrected spectrum of HD 223385 (A3Ia). 
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Fig. 8.— Response curves generated by xtellcor for the various short-wavelength cross- 
dispersed orders of SpeX (8a: order 3; 8b: order 4; 8c: order 5; 8d: order 6; 8e: order 7; 8f. 
order 8). In each plot, the curves shown (from bottom to top) are (i) the observed spectrum 
of an AOV star (HD 223386, V = 6.328), the response curve generated with xtellcor, and 
a theoretical atmospheric transmission spectrum computed for typical conditions on Manna 
Kea, an airmass similar to that of the observations, and a resolving power of 2000. The 
intrinsic stellar H lines of the Paschen and Brackett series are identihed. Small residuals are 
seen in the response curves at the locations of the Pa 11-15 lines. Note that the theoretical 
atmospheric transmission spectrum does not extend below 0.85 pm. 
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Fig. 9.— Final telluric-corrected spectra of an 06.5V((f)) star (HD 206267), generated by 
xtellcor for the various short-wavelength cross-dispersed orders of SpeX (9a: order 3; 96: or¬ 
der 4; 9c: order 5; 9d: order 6; 9e: order 7). In each plot, the spectra are ordered from bottom 
to top as follows: the observed stellar 06.5V((f)) spectrum, the response curve constructed 
by xtellcor from the AOV telluric standard, the telluric-corrected 06.5V((f)) spectrum, and 
the representative theoretical atmospheric transmission spectrum (for reference). 
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Fig. 10.— Final telluric-corrected spectra of a GSIIIa star (HD 16139), generated by 
xtellcor for the various short-wavelength cross-dispersed orders of SpeX (10a: order 3; 106: 
order 4; 10c: order 5; lOd: order 6; lOe: order 7). In each plot, the spectra are ordered from 
bottom to top as follows: the observed stellar GSIIa spectrum, the response curve constructed 
by xtellcor from the AOV telluric standard, the telluric-corrected GSIIIa spectrum, and the 
representative theoretical atmospheric transmission spectrum (for reference). 





























